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NRL  BEAM  PROPAGATION  THEORY  STUDIES  IN  SUPPORT  OF 
SUPERIBEX,  PULSERAD,  RADLAC,  PURE  AND  DELPHI 


OVERVIEW 

This  report  contains  six  short  papers  which  will  appear  in  the 
Proceedings  of  the  Annual  DARPA/SDIO/Services  Charged  Particle  Beam  Review 
which  took  place  at  the  Naval  Postgraduate  School  in  Monterey,  CA  during  18- 
21  September,  1989.  The  papers  describe  electron  beam  propagation  studies 
which  have  been  carried  out  at  NRL  in  support  of  beam  propagation 
experiments  at  several  laboratories. 

Most  of  these  experiments  involve  high  current  propagation  in  air.  In 
some  cases,  the  beam  range  is  extended  by  propagating  the  beam  in  a  reduced- 
der.sity  channel.  SuperlBEX  is  a  new  b  MeV,  10-40  kA  device  located  at  NRL 
which  will  be  used  for  both  beam  stability  and  channel  tracking  studies. 
PULSERAD  is  an  older  1  MeV,  10  kA  diode  at  NRL  which  in  1987  performed  the 
first  successful  channel  tracking  experiments.  RADLAC  is  an  induction 
accelerator  located  at  Sandia  National  Laboratories  which  is  expected  to 
produce  up  to  40  kA  at  an  energy  of  20  MeV.  PURE  is  an  entirely  different 
propagation  concept  which  would  use  an  RF  linac  to  produce  a  train  of 
micron-sized  beam  pulses  at  lower  current  but  very  high  energy.  DELPHI  is 
a  high  altitude  beam  propagation  concept  currently  being  developed  at  Sandia 
which  would  employ  laser-ion  guiding  to  propagate  an  electron  beam  in  a 
diffuse  plasma. 

Some  of  the  papers  included  here  also  describe  predictions  for  the  ATA 
Multi-Pulse  Propagation  Experiment  (MPPE)  which  was  recently  completed  at 
Lawrence  Livermore  National  Laboratory.  NRL  papers  in  the  Annual 
Propagation  Review  which  were  devoted  solely  to  that  experiment  are  included 
in  a  separate  report. 

A  brief  summary  of  each  paper  and  a  list  of  co-authors  for  those  papers 
are  given  below. 

A.  Beam  Propagation  in  Channels:  This  paper  contains  an  overview  of 
NRL  research  on  propagation  in  channels  and  provides  a  brief  summary  of  most 
of  the  papers  described  below  in  this  section  of  the  quarterly  report.  Most 
of  the  work  has  been  in  support  of  the  ATA  Multi-Pulse  Propagation 
Experiment  (MPPE).  Included  are  hose  instability  simulations,  detailed  air 
chemistry  and  channel  physics  calculations  of  ATA/MPPE.  Simulation  studies 
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of  SuperlBEX,  PULSERAD  and  RADLAC  propagation  experiments  are  also 
discussed.  (Hubbard,  Slinker,  Taylor,  Fernsler,  Ali,  Joyce,  Boris) 

B.  Analysis  of  Channel  Tracking  in  SuperlBEX  and  PULSERAD;  SARLAC  was 
used  to  simulate  both  the  1988  PULSERAD  tracking  experiment  and  the  upcoming 
Super-IBEX  experiment.  Expected  trends  were  found;  for  example,  higher  beam 
currents  resulted  in  increased  hose  motion  (shorter  e-folding  distances)  as 
well  as  stronger  tracking  forces  (shorter  tracking  distances).  Other 
interesting  results  emerged  from  the  studies  as  well.  First,  the  general 
features  of  the  PULSERAD  experiments  were  reproduced.  In  particular,  beam 
deflections  and  tracking  distances  were  consistent  with  the  observations. 

The  hose  e-folding  distances  were  comparable  to  the  tracking  distances. 

Also,  the  channel  had  little  apparent  effect  on  the  hose  instability.  While 
there  was  an  apparent  increase  in  hose  growth  in  the  presence  of  a  channel, 
this  effect  was  primarily  due  to  the  reduction  in  scattering;  the  effect 
disappeared  when  the  results  were  properly  scaled  to  the  number  of  betatron 
wavelengths  propagated  by  the  beam.  Second,  Super-IBEX  simulations  showed 
that  beam  deflections  and  tracking  distances  should  be  measurable,  with 
characteristic  hose  distances  being  much  longer  than  tracking  distances. 

That  is,  with  projected  initial  (low  frequency)  perturbations,  the  channel 
may  stabilize  hose.  This  last  effect  is,  in  part,  a  consequence  of  the 
choice  of  the  initial  perturbation  on  the  beam,  which  does  not  couple  to  the 
hose  instability  as  strongly  as  a  high  frequency  perturbation  (PULSERAD). 

The  stabilization  also  may  result  from  higher  return  currents  flowing  at  the 
edge  of  the  channel.  These  simulations  demonstrate  the  importance  of 
limiting  high  frequency  perturbations  in  stability  and  tracking  experiments. 
Similar  predictions  regarding  the  importance  of  initial  beam  perturbation  to 
subsequent  stability  properties  have  been  made  with  regard  to  the  ATA/MPPE 
and  RADLAC.  (Taylor,  Fernsler,  Hubbard,  Slinker) 

C.  Sensitivity  of  Hose  Instability  to  Frequency  of  Initial 

Perturbations  in  Low  and  High  Current  Beams;  Increases  in  the  solenoidal 

guide  field  B  used  in  the  ATA  accelerator  tend  to  reduce  the  growth  of 
z 

high-frequency  BBU  growth  in  the  accelerator  but  enhance  the  generation  of 
low-frequency  sweep  within  the  pulse.  Hose  instability  in  the  propagating 
beam  arises  from  initial  perturbations  generated  by  BBU  and/or  sweep.  To 
examine  the  trade-offs  between  these  two  effects,  a  series  of  SARLAC 
simulations  of  ATA/MPPE  were  performed  with  an  830  MHz  BBU-like  perturbation 
in  the  x-direction  and  a  simple  linear  sweep  in  the  y-direction.  The 
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relative  amplitudes  of  these  perturbations  were  estimated  from  analytical 
models.  In  all  cases,  the  BBU-like  perturbations  grew  more  rapidly  as  the 
beam  propagated  in  air  and  eventually  dominated  even  when  suppressed  by  a  3 
kilogauss  guide  field.  One  encouraging  result  was  that  sweep  amplitudes 
much  larger  than  the  0.01  cm  design  target  for  ATA/MPPE  could  apparently  be 
tolerated.  The  results  strongly  suggest  that  suppression  of  BBU  growth  will 
be  an  important  operational  consideration.  Similar  simulations  were  carried 
out  for  the  higher  current  RADLAC  parameter  regime;  these  simulations  also 
showed  that  high  frequency  perturbations  were  more  dangerous.  One  should  be 
aware  that  high  frequency  perturbations  within  the  beam  lead  to  low 
frequency  oscillations  in  z  because  the  perturbations  couple  to  hose  in  the 
beam  head  where  the  dipole  decay  length  is  short  but  the  betatron  wavelength 
is  long.  (Hubbard,  Slinker,  Taylor) 

D.  Pulse  Decoupling  Using  ATA:  The  channel  created  by  a  multi-pulse 

CPB  burst  becomes  broader  and  shallower  as  the  beam  propagates.  Eventually, 
the  channel  tracking  force  becomes  too  weak  to  guide  the  pulse,  and  the 

geomagnetic  field  deflects  the  beam  out  of  the  channel.  ATA  can  in 

principle  be  used  to  simulate  this  complicated  process  experimentally.  One 
possible  strategy  would  be  to  create  a  channel  with  four  pulses  and  apply  a 
weak  deflecting  magnetic  field  to  the  fifth  pulse.  SARLAC  simulations, 
however,  suggest  that  wall  forces  in  the  ATA  propagation  tank  would  tend  to 
obscure  the  result,  and  the  propagation  distance  required  to  observe  the 
effect  is  relatively  long.  It  may  be  easier  to  study  pulse  decoupling  with 
the  laser-guided  discharge  channels  used  with  SuperlBEX  since  the  tracking 
forces  are  expected  to  be  much  stronger,  the  channel  depth  is  independent  of 
z,  and  the  wall  radius  is  much  larger,  (Fernsler,  Slinker,  Hubbard) 

E.  Electron  Energy  deposition  in  0*.  A  discrete,  time-dependent  energy 

deposition  model  was  used  to  study  high-energy  electron-beam  (100  eV  to  10 
MeV)  deposition  in  0+.  0+  is  expected  to  be  a  major  constituent  in  the  hot 

plasma  channels  created  by  a  PURE-mode  pulse  train.  Secondary  electron 

distributions  were  obtained  by  solving  a  time-dependent,  relativistic 
Boltzmann  equation.  These  distributions  relax  (nonuniformly)  to  steady- 
state  results  from  which  yield  spectra,  production  efficiencies  of  specific 
states,  energy  partitioning,  and  mean  energies  per  electron-ion  pair 
production  (U)  were  computed.  Loss  functions  were  calculated  and  shown  to 
be  in  close  agreement  with  Bethe's  relativistic  equation  for  energies 
greater  than  1  keV.  The  model  predicts  V  is  approximately  72  eV  for  0+  over 
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a  vide  range  of  beam  energies  and  background  ionization  fractions.  The 
effects  of  inner  shell  ionization  and  excitation  were  included  in  the 
deposition  model.  These  effects  result  in  an  increase  in  W  of  approximately 
17%  for  energies  above  10  keV;  this  is  a  direct  consequence  of  the  increased 
ionization  energy.  While  the  present  study,  in  general,  has  focused  on  beam 
sources,  W  was  observed  to  change  for  energies  below  1  keV  if  the  source 
electrons  were  assumed  to  be  completely  stopped  by  the  medium.  (Taylor, 

Ali,  Slinker) 

F.  Transverse  Two-Stream  Instability:  An  intense  relativistic  electron 
beam  propagating  through  a  diffuse  plasma  is  subject  to  a  transverse  two- 
stream  instability.  Analytical  dispersion  relations  have  been  derived 
assuming  an  electrostatic-magnetostatic  approximation  and  spread-mass  beam 
dynamics.  If  a  Bennett-shaped  IFR  channel  is  present,  (the  situation 
expected  for  a  DELPHI  beam),  the  instability  is  hose-like  and  convective, 
and  modest  growth  is  predicted.  However,  if  the  diffuse  plasma  density  is 
constant,  as  expected  for  a  mature  PURE  channel,  the  dispersion  relation 
predicts  absolute  instability  and  rapid  growth.  Simulations  using  the  new 
3-D  ELBA  particle  code  are  consistent  with  the  analytical  model  and  exhibit 
rapid  instability  growth  and  beam  disruption  for  the  case  with  a  constant 
plasma  density.  (Joyce,  Lampe) 
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BEAM  PROPAGATION  IN  CHANNELS* 

R.  F.  Hubbard,  S.  P.  Slinker,  R.  D.  Taylor,  R.  F.  Fernsler 
A.  W.  All,  G.  Joyce  and  P.  Boris 

Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington,  DC  20375 

I.  INTRODUCTION 

Beam  propagation  in  density  channels  is  a  major  focus  of  the  current 
DARPA  experimental  program.  The  ATA  Multi-Pulse  Propagation  Experiment  (MPPE) 
will  attempt  to  demonstrate  stability,  tracking,  and  range  extension  in  the 
channel  formed  by  a  five-pulse  burst  of  10  MeV,  6-8  kA  beams.  Also,  tracking 
experiments  at  NRL  using  the  PULSERAD  and  SuperlBEX  beams  in  laser-guided 
electric  discharge  (LGED)  channels  are  in  progress.  This  paper  provides  an 
overview  of  theoretical  work  at  NRL  in  support  of  these  propagation 
experiments.  More  detail  can  be  found  in  Refs.  1-7. 

II.  ATA  MULTI-PULSE  PROPAGATION  STUDIES 

Overview:  Detailed  predictions  for  an  ATA/MPPE  burst  require  treatment 

of  the  complicated  coupling  between  beam  propagation  and  channel  physics  over 

times  scales  up  to  5  msec.  To  address  these  problems,  we  have  carried  out 

five-pulse  axisymmetric  simulations  which  combine  existing  NRL  propagation  and 

channel  physics  codes. ^  The  results  are  used  to  predict  the  range  of  each 

pulse  and  to  provide  realistic  channels  for  SARLAC  hose  stability  simulations. 

Supporting  these  simulations  are  more  detailed  studies  of  air  chemistry 
2  3 

effects  and  convective  cooling.  We  have  also  studied  the  sensitivity  of 

lead  pulse  hose  instability  growth  to  the  amplitude  of  BBU  and  corkscrew- 

4  5 

induced  perturbations  and  choice  of  emittance  tailoring  method.  We  have 

assessed  the  feasibility  of  an  ATA  pulse-decoupling  experiment.^ 

Typical  parameters  for  these  studies  are  beam  energy  Eq  =  10  MeV,  peak 

current  Iq  =  6  kA,  nominal  beam  radius  a^  =  0.5  cm,  pulse  length  =  33  nsec, 

pulse  separation  t  =  1.25  msec,  and  an  emittance  variation  of  4:1. 

s  1 
Multi-pulse  axisymmetric  propagation  and  channel  dynamics:  In  these 

simulations,  the  SIMMO  particle  code  is  used  for  t  <  to  propagate  the  beam 

and  dump  beam  current  density  Jj}(r,C,z).  Here,  C  is  the  distance  from  the 
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beam  head,  z  is  the  propagation  distance,  and  £,  =  cx  =  ct  -  z.  CHMAIR  II  uses 

the  input  from  SIMMO  to  calculate  beam  and  ohmic  deposition  and  detailed 

air  chemistry  for  x  <  x  <  2x  at  several  z-locations.  HINT  is  then  used  to 

P  P 

calculate  the  long- time-scale  behavior  of  the  channel  (x  <  x  ),  including  the 

s 

effects  of  hydrodynamic  expansion,  thermal  and  convective  cooling,  and 

vibrational  relaxation.  This  generates  a  density  profile,  p(r,z),  which  is 

input  into  SIMMO  for  the  next  pulse  in  the  burst,  and  the  process  is  repeated. 

Six  simulations  were  run  to  determine  the  sensitivity  of  the  channel 

depth  and  the  transported  beam  fluence  to  model  assumptions.  The  fluence, 

Qn(R,z),  is  defined  as  the  transported  beam  charge  for  the  n1^  pulse  within  a 

radius,  R,  of  the  beam  axis  at  a  fixed  location  z.  The  ratio  Q^/Q^  for  R  = 

1.1  cm  and  z  =  6  m  varied  between  1.4  and  1.85,  indicating  modest  range 

extension.  The  channel  depth  at  z  =  0  was  very  sensitive  to  the  assumed 

3  8 

Picone-Boris  form  factor  f  for  convective  cooling;  ’  the  on-axis  density  at 

the  fifth  pulse  was  a  factor  of  two  higher  when  f  was  raised  from  0  to  0.05. 

However, the  predicted  fluence  at  z  =  6  m  changed  by  less  than  15%  because  the 

convectively-cooled  channel  was  significantly  broader.  Changes  in  the  assumed 

chemistry  model  for  SIMMO  and  the  inclusion  of  enhanced  vibrational  cooling 

2 

from  CO^  produced  only  modest  changes  in  the  fluence. 

Chemistry  effects  and  convective  cooling  in  ATA/MPPE;  A  new  chemistry 

model  for  SIMMO  and  SARLAC  was  developed  using  these  same  basic  approach  as  in 

2 

the  standard  "VIPER"  model.  The  new  model  includes  attachment  and  revised 
rate  coeficients  benchmarked  against  detailed  CHMAIR  II  calculations  in  the 
ATA/MPPE  regime.  The  new  mod*1]  gives  similar  axisymmetric  behavior  and 
generates  slightly  more  hose  instability  growth.  The  second  major  focus  of 
the  air  chemistry  studies  was  a  treatment  of  the  transfer  of  the  energy  stored 
in  N£  vibrational  excitations  to  gas  heating.  This  process  occurs  on  the 

millisecond  time  scale  but  can  be  speeded  up  by  adding  a  small  amount  of  C0-. 

8  ^ 
The  Picone-Boris  convective  cooling  form  factor,  which  we  and  others 

have  used  in  MPPE  hydro  simulations  is  a  phenomenological  model  which  has  not 

been  benchmarked  against  full  2-D  hydro  code  results  in  the  appropriate 

parameter  regime.  A  new  2-D  hydro  code  that  treats  the  convective  cooling 

3 

process  directly  has  been  developed.  The  code  runs  on  NRL's  massively- 
parallel  Connection  Machine  and  is  much  faster  than  the  version  used  by  Picone 
and  Boris.  Preliminary  results  suggest  that  the  value  f  =  0.05  used  in  the 
HINT  code  is  reasonable  for  the  level  of  asymmetry  expected  for  ATA/MPPE. 
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Hose  instability  growth  in  MPPE:  SARLAC  simulations^  have  been  carried 
out  in  the  channels  generated  by  the  SIMMO-CHMAIR-BINT  simulation  sequence 
described  above.  Hose  amplitudes  in  the  beam  body  grow  om  an  assumed 
initial  level  of  0.01  cm  to  0.5-0. 7  cm  in  6  m  of  propagation.  Similar  hose 
amplitudes  were  observed  in  the  absence  of  a  channel. 

Increasing  the  guide  field  B  in  ATA  increases  low  frequency  sweep 

z  4 

displacements  but  suppresses  high  frequency  BBU-induced  perturbations.  The 
effect  of  this  tradeoff  on  hose  instability  was  investigated  using  SARLAC  by 
initiating  a  BBU-like  mode  in  x  and  a  linear  sweep  in  y.  The  BBU-induced  mode 
was  much  more  unstable  in  the  simulations,  while  relatively  large  sweep 
amplitudes  could  be  tolerated.  In  one  set  of  SARLAC  runs,  increasing  Bz  from 
1  to  3  kG  caused  hose  amplitudes  at  z  =  6  m  to  drop  by  almost  a  factor  of  3. 

Other  SARLAC  simulations  considered  the  effects  of  different  possible 
emittance  tailoring  techniques.^  Beam  radius  and  emittance  profiles  were 
generated  by  FRIEZR  for  various  tailoring  schemes,  and  the  results  were  used 
to  initialize  SARLAC.  A  multi-foil  tailoring  cell  simulation  produced  hose 
amplitudes  in  excess  of  1  cm  in  4.8  m  of  propagation,  while  a  similar  run 
using  a  5  mtorr  IFR  cell  grew  to  only  0.1  cm. 

Pulse  decoupling  experiment  for  ATA:6  SARLAC  was  used  to  investigate  the 
feasibility  of  studying  pulse  decoupling  experimentally  on  ATA  by  applying  a 
weak  deflecting  guide  field  just  before  the  last  pulse.  The  simulation 
results  suggest  that  the  phenomenon  would  likely  be  obscured  by  wall  forces. 

Summary:  The  ATA/MPPE  beams  are  likely  to  propagate  with  only  moderate 
hose  growth  if  the  stringent  beam  conditioning  goals  are  met.  Range  extension 
effects  should  be  modest  but  observable. 

III.  SUPERIBEX  AND  PULSERAD  TRACKING  STUDIES 

PULSERAD  stability  and  tracking  simulations: ^  A  data  base  now  exists  for 
the  1988  NRL  tracking  experiments  which  used  the  1  MeV  PULSERAD  beam.  SARLAC 
simulations  were  performed  with  a  weakly-tailored  beam  with  the  estimated 
experimental  parameters.  The  hose  amplitudes  were  somewhat  higher  in  a 
centered  density  channel  than  in  full  density  air,  but  the  difference  could  be 
attributed  entirely  to  scattering  effects.  Increasing  Iq  appeared  to 
destabilize  the  beam.  The  tracking  distance  was  typically  30  cm.  All  of 
these  results  are  consistent  with  the  experimental  data. 
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SuperlBEX  s lability  and  tracking  simulations:  SARLAC  was  used^  to  model 
the  4.5  MeV  SuperlBEX  tracking  experiments  currently  in  progress.  The 
simulations  assumed  moderate  emittance  tailoring  and  relatively  large 
amplitude,  but  low  frequency  initial  hose  perturbations.  Peak  simulation  beam 
currents  varied  between  10  and  40  kA.  Hose  amplitudes  increased  with  Iq  and 
were  smaller  in  the  presence  of  a  channel  than  in  full  density  air,  in 
contrast  to  the  PULSERAD  result.  Tracking  distances  were  predicted  to  be  only 
20-30  cm  in  the  presence  of  an  offset  channel. 

A  separate  series  of  SARLAC  simulations  modeled  propagation  in  a  10  m 
long  tank  in  uniform  0.5  atm  air  and  in  a  centered,  0.1  atm  density  channel. 
The  beam  was  well-conditioned,  as  might  be  produced  by  a  two-stage  IFR/active- 
wire  conditioning  cell  currently  being  considered  for  future  experiments.  The 
perturbations  were  the  same  as  used  in  the  RADLAC  simulations  described  in 
Ref.  4.  Hose  amplitudes  grew  to  only  0.2-0. 3  cm  in  the  SuperlBEX  simulations. 

Summary:  PULSERAD  simulation  results  were  consistent  with  the  data  from 
the  successful  1988  tracking  experiment.  Simulations  of  the  current  SuperlBEX 
tracking  experiment  predict  short  tracking  distances  and  moderately-strong 
hose  growth.  Future  SuperlBEX  experiments  with  well-conditioned  beams  could 
exhibit  "stable"  propagation  to  ranges  beyond  10  m. 
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Analysis  of  Channel  Tracking  in  Super-IBEX  -and  PULSKRAD* 

R.D.  Taylor,  R.F.  Fernsler,  R.F.  Hubbard,  and  S.P.  Slinker 
Beam  Physics  Branch,  Plasma  Physics  Division 
Naval  Research  Laboratory,  Washington,  D.C.  20375 

SARLAC  is  used  to  simulate  both  the  1988  PULSERAD  tracking  experiment 
and  the  upcoming  Super-IBEX  experiment.  From  a  theoretical  (and  simulation) 
perspective,  key  questions  about  these  experiments  include:  Do  the  simulations 
predict  a  tracking  force  and,  if  so,  on  a  scale  that  is  observable?  Does  the 
existence  of  a  channel  affect  the  stability  properties  of  the  propagating 
beam?  Several  interesting  results  emerge  from  these  studies.  First,  the 
general  features  of  the  PULSERAD  experiments  are  reproduced.  In  particular, 
beam  deflections  and  tracking  distances  are  consistent  with  the  observations. 
Also,  the  channel  has  little  apparent  effect,  on  the  hose  instability.  Second, 
Super-IBEX  simulations  show  that  beam  deflections  and  tracking  distances 
should  be  measurable.  In  addition,  with  projected  initial  (low  frequency) 
perturbations,  the  channel  may  stabilize  hose.  This  last  effect  is,  in  part,  a 
consequence  of  the  choice  of  the  initial  perturbation  on  the  beam.  Similar 
predictions  regarding  the  importance  of  initial  beam  perturbation  to 
subsequent  stability  properties  have  been  made  elsewhere.* 

The  PULSERAD  simulations  assume  a  weakly  tailored  (1.2: 1.0  over  8  ns), 

1  MeV  beam  having:  1)  a  4G  ns  triangular  temporal  profile  with  a  20  ns 
risetime  and  20  ns  falltime;  2)  a  nominal  radius  of  2.0  cm;  3)  a  variable 

frequency  (420  MHz  at  the  head  to  102  MHz  at  the  tail)  initial  hose 

perturbation  with  a  2  mm  amplitude;  and  4)  currents  of  7  and  10  kA.  The 
channel  characteristics  are:  1)  a  minimum  density  equal  to  0.1  times  that  of 

ambient  air  (pQ);  2)  a  4.0  cm  radius;  and  3)  an  offset,  where  applicable,  of 

1.4  cm.  The  beam  is  injected  into  full  density  air  (with  the  scattering  turned 
on  and  off  to  delineate  stability  properties),  an  on-axis  channel,  and  an 
off-axis  channel. 

The  Super-IBEX  simulations  assume  a  moderately  tailored  (2. 5:1.0  over 
7  ns),  5  MeV  beam  having:  1)  a  35  ns  temporal  profile  with  a  15  ns  linear 

*  Work  supported  by  the  Defense  Advanced  Research  Projects  Agency,  ARPA  Order 
No.  4395,  Amendment  No.  80,  monitored  by  the  Naval  Surface  Warfare  Center. 

**  Address:  Berkeley  Research  Associates,  Springfield,  VA  22150 
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rise,  5  ns  plateau,  and  15  ns  fall;  2)  a  nominal  radius  of  1.5  cm;  3)  a 
constant.,  low  frequency  (33  MHz)  initial  hose  perturbation  with  a  2  mm 
amplitude;  and  4)  currents  of  10,  20,  30,  and  40  kA.  The  channel 
characteristics  are:  1)  a  minimum  density  equal  to  0.2pq;  2)  a  4.0  cm  radius; 
and  3)  an  offset,  where  applicable,  of  1.5  cm.  The  beam  was  injected  into  full 
density  air,  reduced  density  (0.2pQ)  air,  an  on-axis  channel,  and  offset  from 
the  channel  axis. 

In  both  cases,  tracking  distances  and  hose  e-folding  distances  are 
computed.  The  channel/beam  offset  is  initiated  in  the  y-direction  so  hose 
e-folding  distances  are  obtained  from  the  x-centroid  motion.  Sample  results 
for  the  20  ns  slice  of  each  beam  (spanning  the  designated  currents)  are  shown 
in  the  following  table: 


PULSERAD 

Propagation 

Medium 

Tracking 

Distance  (cm) 

Hose  e-folding 
Distance  (cm) 

Full  air 

— 

29  -  34 

Centered  channel 

— 

23  -  27 

Offset  channel 

28  -  33 

18  -  27 

Super-IBEX 

Propagation 

Medium 

Tracking 

Distance  (cm) 

Hose  e-folding 
Distance  (cm) 

Full  air 

— 

40  -  96 

0.2  atm  air 

— 

40  -  80 

Centered  channel 

— 

>  200 

Offset  beam 

21  -  28 

144  -  164 

The  simulation  results  follow  the  generally  expected  trends.  For  example, 
higher  beam  currents  result  in  increased  hose  motion  (shorter  e-folding 
distances)  as  well  as  stronger  tracking  forces  (shorter  tracking  distances). 

The  PULSERAD  simulations  show  tracking  distances  comparable  to  hose  e- 
folding  distances.  The  centered-channel  cases  appeared  to  be  somewhat  more 
unstable  than  the  full-density-air  cases.  Both  results  are  consistent  with 
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Fig.l.  7  kA  PULSZRAD  beam  injected  into 
full  density  air;  y-centroid 
(dashed  line),  x-centroid  (solid 
line). 


Fig. 2.  7  ka  PULSZRAD  beam  injected 
into  offset  channel;  same  as 
in  Fig.  1,  channel  edges  and 
axis  (dot-dashed  line). 


4 - - 
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Fig. 3.  30  kA  Super-IBZX  beam  injected 
along  channel  axis;  same  as  in 
Fig.  2. 


Fig. 4.  30  kA  Super-IBZX  beam 

injected  at  an  offset  from 
channel;  same  as  in  Fig.  2. 


2 

those  found  in  the  experiment.  A  comparison  of  Figs.  1  and  2  shows  the  effect 
of  an  offset  channel  on  the  beam.  Not  shown  in  the  figures  are  cases  for  the 
beam  in  a  centered  channel  and  the  beam  in  air  with  the  scattering  turned  off. 
These  latter  cases  have  nearly  identical  current  and  radius  profiles  and 
similar  hose  instability  amplitudes.  Thus,  the  apparent  increase  in  hose 
growth  in  the  presence  of  the  channel  is  primarily  due  to  the  reduction  in 
scattering.  When  hose  growth  is  properly  scaled  to  the  number  of  betatron 
wavelengths  propagated  by  the  beam,  the  channel  has  virtually  no  effect.  Over 
90  cm,  the  beam  deflection  for  these  cases  is  comparable  to  the  x-deflection 
seen  in  Fig.  2  which,  in  turn,  is  greater  than  that  shown  in  Fig.  1. 

The  Super-IBEX  simulations  show  tracking  distances  much  shorter  than  the 
characteristic  hose  distances.  Figure  4  is  a  dramatic  example  of  the  effect  of 
the  tracking  force.  A  beam  injected  off  of  the  channel  axis  is  drawn  toward 
the  channel,  overshoots,  rides  the  far  edge,  is  drawn  back  into  the  channel, 
and  begins  to  oscillate.  A  beam  injected  along  the  channel  axis  is  stabilized 
over  a  substantial  distance  (225  cm);  see  Fig.  3.  As  seen  in  the  table,  the 
beam  is  much  less  unstable  in  the  centered  channel  than  in  full  density  air 
even  though  it  has  traveled  more  betatron  wavelengths.  There  are  two  reasons 
for  this  apparent  stabilization.  First,  the  higher  currents  in  Super-IBEX 
result  in  much  more  return  current  flow  at  the  edge  of  the  channel  which  is  a 
well-known  stabilizing  effect.  In  addition,  the  Super-IBEX  runs  assume  a 
constant,  low  frequency  perturbation,  representative  of  that  expected  in  the 
experiment.  These  perturbations  do  not  couple  to  the  hose  instability  as  well 
as  the  higher  frequency  modes  used  to  initialize  the  PULSERAD  simulations. 
These  simulations  are  a  clear  demonstration  of  the  importance  of  limiting  high 
frequency  perturbations  in  stability  and  tracking  experiments.^- 


1.  R.  Hubbard,  S.  Slinker,  and  R.  Taylor,  "Sensitivity  of  Hose  Instability  to 
Frequency  of  Initial  Perturbations  in  Low  and  High  Current  Beams",  these 
proceedings. 

2.  D.  Murphy,  "Evidence  of  Hose  Instability  in  Density  Channels",  these 
proceedings. 

3.  R.  Meger,  J.  Antoniades,  and  T.  Peyser,  private  communication. 
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SENSITIVITY  OF  HOSE  INSTABILITY  TO  FREQUENCY  OF  INITIAL  PERTURBATION 

IN  LOW  AND  HIGH  CURRENT  BEAMS* 

R.  F.  Hubbard,  S.  P.  Slinker  and  R.  D.  Taylor** 

Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington,  DC  20375 

I.  INTRODUCTION 

Most  studies  of  resistive  hose  instability  in  propagating  electron  beams 

have  assumed  that  the  frequency  spectrum  of  the  initial  perturbation  cannot  be 

experimentally  controlled.  However,  Fawley  has  pointed  out  that  for  the  ATA 

Multi-Pulse  Propagation  Experiment  (MPPE),  the  relative  amplitudes  of  low 

frequency  sweep  and  high  frequency  (BBU-induced)  perturbations  are  sensitive 

to  the  guide  field,  B  ,  in  the  accelerator.  Increasing  B  reduces  the  growth 

z  z 

of  BBU  but  enhances  corkscrew-induced  sweep. 

We  have  used  the  SARLAC  hose  instability  code^  to  determine  which  class 

of  perturbations  is  the  more  dangerous.  SARLAC  models  beam  propagation  in  the 

atmosphere  using  the  doppler-shif ted  coordinates  z  and  C  =  ct  -  z.  Here,  z  is 

the  propagation  distance  in  air  and  C  is  the  distance  from  the  beam  head.  A 

BBU-like  830  MHz  perturbation,  X(C,z=0),  is  imposed  in  the  x-direction  using 

2 

an  asymptotic  BBU  growth  model  from  Caporaso.  A  low  frequency  (LF)  sweep 

perturbation  is  imposed  in  the  y-direction.  In  general,  the  hose  growth  in 

the  x-direction  dominates  even  when  its  initial  amplitude  is  substantially 

lower,  suggesting  that  ATA  should  be  operated  with  a  high  B  .  Similar  results 

3  z 

have  been  obtained  by  Feinstein  and  Keeley.  HF  perturbations  also  generate 
more  instability  growth  in  RADLAC  simulations.  These  HF  perturbations 
initiate  hose  growth  in  the  "neck"  of  the  beam  where  the  dipole  decay  length 
is  relatively  short.  Since  the  local  betatron  wavelength  is  usually  long  in 
this  region,  an  HF  perturbation  may  generate  low  frequency  oscillations  in  z. 

II.  ATA  HOSE  SIMULATIONS 

2 

BBU  perturbation  model;  Following  Caporaso,  the  beam  exits  the 
accelerator  with  a  displacement  X(C,z=0)  =  Xpexp(KNI^(  C)Z^(Op/Bz)sin(&)0C/c) . 
Here,  I.  is  the  beam  current  in  kA,  the  cavity  impedance  Z.  =  30  5,  <of./2n  = 
8.3x10  sec  ,  the  number  of  cavities  N  =  50,  B  =  1-3  kG,  Xq  =  10  cm  and 
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the  coefficient  K  =  1.16x10  kG-sec(kA-ohm)  .  The  beam  is  assumed  to  have 
an  energy  of  10  MeV  with  the  current  ramping  to  its  nominal  value  Iq  =  6  or  8 
kA  over  a  distance  Cr  =  360  cm.  The  injected  radius  and  emittance  profiles 
follow  the  form  shown  in  Fig.  4  of  Ref.  4  which  are  generated  from  a  FRIEZR 
simulation  of  a  5  mtorr  passive  IFR  conditioning  cell. 

Corkscrew  or  sweep  perturbation  model;  Low  frequency  sweep  or  corkscrew 
arises  from  the  coupling  between  energy  variations  Ay  within  the  pulse  and 
field  or  alignment  errors.  This  produces  a  phase  advance,  &4>  =  (Ay/y) J“kcdz, 
where  the  cyclotron  wavenumber  k  is  proportional  to  B  .  For  <<  1,  the 

C  2 

phase  advance  and  sweep  amplitudes  are  proportional  to  Bz-  In  SARLAC,  this 
effect  is  modeled  by  imposing  an  initial  perturbation  in  the  y-direction  given 
by  Y(C,z=0)  =  aB^C*  The  coefficient  a  is  chosen  to  be  3xl0~^  (kG)-*. 

Results:  Four  long  SARLAC  simulations  were  run  with  the  parameters 
described  above.  Cases  Al,  A2  and  A4  were  for  Iq  =  6  kA  and  Bz  =  1,  1.5  and 
3  kG,  respectively,  while  Case  A3  used  an  8  kA  beam  with  a  1.5  kG  guide  field. 
Results  are  summarized  in  the  table  below.  The  displacements  X  and  Y  are 
taken  in  the  beam  tail  at  C  =  750  cm  (25  nsec)  and  are  tabulated  at  injection 
(z  =  0)  and  at  z  =  5.4  m. 


Case 

*o 

Bz 

X(0) 

Y(0) 

X(5.4m) 

Y(5.4m) 

Al 

6  kA 

1.0  kG 

0.010  cm 

0.022  cm 

0.31  cm 

0.12  cm 

A2 

6 

1.5 

0.0035 

0.034 

0.20 

0.12 

A3 

8 

1.5 

0.010 

0.034 

0.67 

0.14 

A4 

6 

3.0 

0.0006 

0.068 

0.11 

0.12 

Figure  1 

plots 

the  initial 

perturbations  X(£,0) 

and  Y(C,0)  for  Case  A2, 

showing  that 

the  HF 

perturbation 

(solid  curve)  is  an 

order  of 

magnitude  lower 

than  the  sweep  perturbation  (dashed  curve).  However,  Fig.  2  shows  that  at  z  = 
6  m,  the  HF  modes  induced  in  the  x-direction  have  surpassed  those  in  the  y- 
direction.  Increasing  Iq  to  8  kA  results  in  a  substantial  increase  in  X,  both 
because  the  BBU-induced  initial  perturbation  is  larger  and  because  the  higher 
current  beam  travels  more  betatron  wavelengths.  The  resulting  hose  amplitudes 
are  shown  in  Fig.  3.  The  results  suggest  that  ATA  should  be  operated  with  a 
relatively  high  guide  field  to  suppress  BBU  and  that  relatively  large  sweep 
amplitudes  may  be  tolerated  without  initiating  serious  instability  growth. 
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Fig  1.  Initial  dis¬ 
placement  X(C,0)  and 
Y(C*0)  (solid  and 
dashed  lines):  Case  A2 


Fig  2.  X(0  and  Y(t) 
at  z  =  6m  for  Case  A2. 
BBU-induced  hose  is 
stronger  (solid  line). 


(  (an) 


Fig  3.  X(C)  and  Y(C)  at 
z  =  8m  or  Case  A3.  1q  is 
raised  to  8  kA;  beam  is 
more  unstable. 


III.  RADLAC  HOSE  INSTABILITY  SIMULATIONS 

Overview:  Both  HF  and  LF  perturbations  are  produced  in  RADLAC  even 

though  BBU  is  thought  to  be  unimportant.  A  series  of  SARLAC  simulations  were 

performed  using  a  0.02  cm,  830  MHz  HF  perturbation  in  x  and  a  0.2  cm  (at 

C  =  900  cm)  sweep  perturbation  in  y.  Nominal  beam  parameters  in  the 

simulations  were  Iq  =  25  kA,  Tq  =  41  (ramped  in  some  cases),  beam  radius 

a^  *  2  cm,  rise  length  =  300  cm  and  normalized  emittance  taper  Et  =  2-4. 

RADLAC  simulation  results:  Six  simulations  were  performed  as  described 

in  the  table  below.  E  .  and  E  define  the  range  of  the  energy  ramp,  and 

min  max 

X  and  Y  are  the  maximum  hose  amplitudes  at  C  =  900  cm  and  z  <  12  m. 
max  max 


Case 

E  . 
min 

E 

max 

Comment 

X 

max 

Y 

max 

R1 

20  MeV 

20  MeV 

4 

0.5  cm 

0.5  cm 

R2 

20 

20 

4  a,=1.5cm 

2.5 

0.8 

R3 

5 

20 

2  b  - 

>  20 

>  5 

R4 

5 

20 

4 

4.5 

1.5 

R5 

5 

20 

4  YQ=0.02cm 

4.5 

0.6 

R6 

10 

20 

4  Faster  y-ramp 

1.1 

0.4 

As  in  the  ATA  simulations,  the  high  frequency  modes  dominate  even  though 
they  are  initiated  at  a  smaller  amplitude.  Comparing  Cases  R1  and  R3,  it  is 
apparent  that  relying  on  the  natural  tailoring  which  comes  about  from  the 
energy  ramp  may  lead  to  unacceptably  large  hose  amplitudes.  This  is,  in  part, 


15 


*>.  T»  (cm) 


because  the  head  is  so  hot  that  it  is  quickly  lost,  Leaving  behind  a  poorly- 
tailored  beam.  Even  when  the  4:1  emittance  variation  is  restored  to  a  beam 
with  a  Y-ramp  (Case  R4),  the  beam  is  more  unstable  than  in  the  constant  energy 
case.  Figure  4  plots  X(z)  and  Y(z)  at  (  =  675  cm  for  Case  R4,  showing  an 
initial  damping  of  the  LF  mode  followed  by  an  eventual  coupling  to  the  faster¬ 
growing  HF  mode.  Hose  amplitudes  versus  C  at  z  =  120  cm  are  shown  for  this 


case  (Fig.  5)  and  for  the  more  unstable  weakly-tapered  Case  R3  (Fig.  6). 


t  (on)  {  (cm)  {  (cm) 


Fig  4  X(z)  and  Y(z)  at  Fig  5.  X(C)  and  Y(C)  at  Fig.  6.  X(C)  and  Y(C)  at 

C  =  675cm  for  RADLAC  z  =  12  m  for  Case  R4.  z  =  12  m  for  Case  R4 

Case  R4.  Note  initial  HF  mode  (solid  line)  (weaker  taper).  Note 

decay  in  Y  (dash  line).  dominates.  shorter  pulse  length. 


IV.  CONCLUSIONS  AND  REFERENCES 


For  both  ATA  and  RADLAC,  high  frequency  perturbations  appear  to  couple 
more  strongly  to  the  resistive  hose  instability  and  should  be  suppressed  if 
possible  even  at  the  expense  of  higher  sweep  amplitudes. 
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PULSE  DECOUPLING  USING  ATA*  . 


R.  Fernsler,  S.  Slinker,  and  R.  Hubbard 
Beam  Physics  Branch,  Plasma  Physics  Division 
Naval  Research  Laboratory,  Washington,  D.C. 


INTRODUCTION 

As  the  lead  pulse  propagates,  its  energy  and  geomagnetic  gyroradius 
decrease  while  its  radius  increases.  The  density  channel  formed  by  the  pulse 
thus  becomes  progressively  curved  but  broader  and  more  shallow.  As  a 
consequence,  a  follower  beam  pulse  is  predicted^  to  decouple  from  the  lead- 
pulse  channel  as  the  channel  tracking  force  becomes  weak  (in  the  beam  head) 
relative  to  the  geomagnetic  force.  The  follower  pulse  then  becomes  the  new 
lead  pulse  and  constructs  its  own  channel.  Although  predictable  pulse 
decoupling  is  a  critical  applications  issue,  it  has  received  only  cursory 
analytic  attention  and  no  experimental  or  numerical  verification. 

In  this  paper  we  consider  using  ATA  for  studying  pulse  decoupling.  A 
series  of  four,  closely  spaced  beam  pulses  would  be  used  to  create  a  density 
channel  in  atmospheric  air.  Prior  to  the  fifth  pulse,  a  fixed  external 
magnetic  field  -  1  Gauss  would  be  applied.  Because  the  density  channel 
weakens  with  range  z,  the  tracking  force  would  weaken  until  the  fifth  pulse 
ultimately  left  the  channel  created  by  its  predecessors. 

ISSUES 

Pulse  decoupling  is  a  complicated  process  involving  channel  evolution  and 

the  mutual  interaction  of  three  forces:  an  external  (geomagnetic/centrifugal) 

deflection  force  generated  by  the  geomagnetic  field,  channel  curvature,  and 

2 

pulse-to  pulse  energy  differences;  a  channel  tracking  force  arising  from 
beam-generated  ionization,  air  chemistry,  and  the  depression  in  channel  gas 
density;  and  a  body  coupling  force  that  causes  the  beam  body  to  follow  the 
head  (and  also  drives  the  resistive  hose  instability).  Existing  codes  should 
be  capable  of  modeling  all  of  these  forces  with  only  minor  modifications. 
Experimental  verification,  however,  is  more  difficult. 

An  ideal  pulse-decoupling  experiment  should  address  several  key  issues: 
(1)  How  strong  can  the  external  force  be  before  channel  tracking  and  pulse 

♦Work  supported  by  the  Defense  Advanced  Research  Projects  Agency,  ARPA  Order 
No.  4395,  Amendment  80,  monitored  by  the  Naval  Surface  Warfare  Center 
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coupling  are  destroyed?  (2)  What  are  the  effects  of  .beam  tilt  on  the  hose 
instability  and  hole-boring?  (3)  Hov  much  emittance  tailoring,  for  hose 
suppression,  can  be  tolerated?  (4)  Is  pulse  decoupling  a  smooth  and 
predictable  process,  and  when  and  where  does  it  occur?  And  (5),  to  what 
extent  does  an  external  force  enhance  nose  erosion?  Answering  any  of  these 
questions  requires  a  beam  that  remains  stable  and  tightly  pinched  as  it  comes 
to  equilibrium  with  the  various  deflection  forces. 

SIMULATION 

To  assess  the  potential  of  ATA  for  studying  pulse  decoupling,  we  used  the 
SARLAC  code  to  simulate  a  single  ATA  pulse  propagating  into  a  preexisting 
channel.  The  beam  parameters  at  injection  were:  peak  current  6  kA;  half¬ 
current  radius  R^^  =  cm  *n  t^le  beam  body;  emittance  tailoring  of  4-to-l 
over  325  cm;  matching  current  I  =  3  kA;  and  displacement  x  =  0.75  cm  from 
the  channel  axis.  The  density  channel  was  Gaussian  with  a  radius  of  1.5  cm 
and  on-axis  density  of  0.2  atm.  The  channel  was  uniform  in  z  and  centered 
within  a  conducting  pipe  of  radius  b  *  32  cm.  A  fixed  deflection  force  of  3 
Gauss  was  applied.  The  choice  of  a  large  pipe  radius  and  large  deflection 
force  is  explained  below. 

RESULTS 

Figure  1  shows  the  beam  centroid  (x,y)  as  a  function  of  distance  behind 
the  beam  head,  C  =  ct-z,  at  several  propagation  distances  z  in  the  absence  of 
a  channel.  The  beam  tilts  in  response  to  the  wall  forces  which  become  strong 
in  the  body  and  retard  its  off-axis  motion.  The  body  consequently  remains 
closer  to  the  pipe  axis  than  does  the  head.  [The  apparently  modest  drift  of 
beam  slices  at  t  <  50  cm  is  an  artifact  of  beam  expansion  and  scrape-off. J  A 
significant  reduction  in  the  pipe  radius  or  deflection  force  is  undesirable 
because  the  wall  forces  would  then  obscure  the  deflection  force.  Increasing 
the  pipe  radius  to  b  >  1  m  would  virtually  eliminate  the  wall  forces. 

Figure  2  plots  the  beam  centroid  (x,y)  with  the  channel  present,  and 
illustrates  the  effect  of  the  channel  tracking  force.  This  force,  together 
with  the  wall  force,  exceeds  the  external  force  of  3  Gauss  at  C  >  120  cm,  and 
pulls  the  beam  body  into  the  channel.  The  beam  head,  however,  is  pulled  out 
of  the  channel  by  the  external  deflection  force.  Observe  that  the  beam  tilt 
is  increasing  at  z  =  3  m,  indicating  that  the  beam  has  not  yet  reached 
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Fig.  1.  Beam  centroid  (x,y)  at  z  -  0,  1.5,  2.25,  3  m  without  a  channel. 


equilibrium.  Moreover,  whether  the  beam  body  remains  in  the  channel  or 
follows  the  ejected  head  is  unclear.  Hence,  this  simulation  demonstrated 
channel  tracking  but  not  pulse  decoupling.  Range  extension  is  evident  in 
Fig.  3  in  that  the  beam  expands  less  when  propagating  in  the  channel. 

CONCLUSION 

Simulation  of  a  simple  pulse-decoupling  experiment  for  ATA  revealed 
several  difficulties,  the  principal  ones  being  the  wall  forces  and  the  need 
for  the  beam  to  propagate  a  long  distance  before  pulse  decoupling  can  be 
observed.  A  potential  solution  to  the  latter  problem  is  to  inject  the  beam 
on  axis  in  the  channel  and  turn  on  the  external  deflecting  force  gradually  in 
z.  The  beam  should  then  remain  in  near  equilibrium  with  the  external, 
tracking,  and  coupling  forces. 
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Fig.  3.  ri/2  (a)  and  with  (b)  a  density  channel.  Plots  at  z  *  0, 

0.75,  1.5,  2.25,  and  3  m. 
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Electron  Energy  Deposition  in  0+  * 

R.D.  Taylor,  A.W.  Ali,  and  S.P.  Slinker 
Beam  Physics  Branch,  Plasma  Physics  Division 
Naval  Research  Laboratory,  Washington,  D.C.  20375 

A  discrete,  time-dependent  energy  deposition  model  is  used  to  study 
high-energy  electron-beam  (100  eV  to  10  MeV)  deposition  in  0+.  Details  of  the 
deposition  code  and  applications  to  other  gases  have  been  reported 
previously^  ^  and  published  elsewhere. ^  ^  Secondary  electron  distributions  are 
obtained  by  solving  a  time-dependent,  relativistic  Boltzmann  equation.  These 
distributions  relax  to  steady-state  results  from  which  yield  spectra, 
production  efficiencies  of  specific  states,  energy  partitioning,  and  mean 
energies  per  electron-ion  pair  production  (V)  are  computed.  Loss  functions  are 
calculated  and  compared  to  Bethe's  relativistic  equation.^  The  model  predicts 
W  is  approximately  72  eV  for  0+  over  a  wide  range  of  beam  energies  and 
background  ionization  fractions.  The  effects  of  inner  shell  ionization  and 
excitation  have  been  included  in  the  deposition  model.  These  effects  result  in 
an  increase  in  V  of  approximately  17%  for  energies  above  10  keV.  While  the 
present  study  has  focused  on  beam  sources,  W  is  observed  to  change  for 
energies  below  1  keV  if  the  source  electrons  are  assumed  to  be  completely 
stopped  by  the  medium. 

The  Boltzmann  equation  describes  a  spatially  homogeneous  electron  beam 
impinging  upon  a  single  component  gas.  Electric  and  magnetic  field  effects  due 
to  the  beam  are  neglected.  The  effects  of  inelastic  and  ionizing  collisions 
are  accounted  for,  however,  elastic  and  superelastic  collisions  between  the 
electrons  and  background  gas  are  neglected.  The  gas  density  and  source  term 
are  time  independent.  Energy  loss  to  the  plasma  electrons  is  included  by  means 
of  a  loss  funtion.  Coupling  to  the  radiation  field  is  neglected. 

For  deposition  in  0+,  energy  loss  to  the  2p  ZD°  and  2pJ  ^P°  states,  all 

3  4*  4- 

n  =  3  states,  2  effective  Rydberg  series  (n  >  4)  states,  and  the  3P  0**  state 
is  considered.  Analytic  forms  for  the  relevant  excitation  and  ionization  cross 
sections  have  been  obtained.  A  detailed  discussion  of  these  forms,  their 
reduction  to  well-known  theoretical  results,  and  their  relative  accuracy  in 

*  Work  supported  by  the  Defense  Advanced  Research  Projects  Agency,  ARPA  Order 
No.  4395,  Amendment  No.  80,  monitored  by  the  Naval  Surface  Warfare  Center. 

**  Address:  Berkeley  Research  Associates,  Springfield,  VA  22150 
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comparison  to  existing  measurements  and  calculations  is  given  elsewhere.  ’ 

The  K-shell  ionization  threshold  is  565  eV  and  it  is  assumed  that  all  Auger 
electrons  are  produced  at  475  eV. 

18  -2  -1 

Results  are  generated  for  a  beam  flux  of  1.99x10  cm  sec  ,  a 

19  -3 

background  density  of  2.46x10  cm  ,  and  various  fractional  ionizations  (zero 
unless  otherwise  noted).  Figure  1  shows  the  loss  function,  L(T),  (neglecting 
inner-shell  effects)  and  its  components  for  energies  up  to  10  MeV  for  0+.  Most 
of  the  energy  goes  into  producing  secondary  electrons,  while  loss  to 
excitation  is  significant  only  below  the  ionization  threshold.  Also  shown  is  a 
comparison  with  Bethe's  relativistic  loss  function,  Lg(T).  For  energies 
greater  than  1  keV,  L(T)  and  Lg(T)  are  in  close  agreement.  In  principle, 
energy  loss  of  an  electron  traversing  a  material  is  lessened  because  of 
polarization  of  the  medium.  This  density  effect  is  small  for  the  energies  in 
the  present  studies.  The  inclusion  of  inner-shell  effects  and  the  resultant 
changes  to  the  loss  function  are  shown  in  Fig.  2.  For  example,  for  energies 
greater  than  1  keV,  the  inner  shell  contribution  exceeds  loss  to  excitation. 

Average  excitation,  ionization,  and  secondary  energies  are  presented  in 
Fig.  3.  There  is  very  little  change  in  I  and  Eg  for  energies  greater  than  ~ 

100  eV.  Asymptotically,  I  =  35.1  eV  and  =  20.6  eV  for  0+.  K-shell 
ionization  increases  I,  also  shown  in  Fig.  3. 

The  steady-state  values  of  V  (with  and  without  inner-shell  ionization) 
for  beam  energies,  Tfa,  ranging  from  100  eV  to  10  MeV  are  shown  in  the 
following  table: 


Tb  (ev) 

102 

103 

io4 

105 

106 

io7 

m 

80.5 

72.3 

72.3 

72.2 

71.3 

70.0 

V(eV) :  L  &  K 

— 

74.2 

82.0 

84.2 

84.9 

84.9 

These  values  are  nearly  constant  over  an  large  energy  range.  The  near 
constancy  of  V  at  energies  >  1  keV  is  well-known.  The  increase  in  V  which 
accompanies  the  addition  of  K-shell  ionization  is  a  direct  conequence  of  the 
increased  I.  Figure  4  shows  a  comparison  between  V  for  beam  electrons  and 
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AVERAGE  EHEI1GY 


for  completely-stopped  electrons.  There  is  very  little  difference  for  energies 

greater  than  10  keV,  but  significant  differences  below. 

Other  results  obtained  from  the  present  analysis  and  discussed  in  detail 
9 

elsewhere  include  the  following: 

Relaxation  of  the  distribution  function  to  a  steady-state  is  non- 

uniform,  i.e.  energy-dependent.  Intermediate  energies  (-  100  eV)  relax  first, 

followed  by  the  lower  part  of  the  spectrum  (-  10  eV,  but  greater  than  the 

lowest  excitation  energy),  and,  finally,  the  high  energies  (>  1  keV). 

Production  efficiencies  (per  electron-ion  pair  created)  for  excited 

states  are  nearly  constant  for  the  beam  energies  considered. 

The  results  show  some  sensitivity  to  background  ionization  fraction.  In 

particular,  the  distribution  function  is  insensitive  to  changes  in  the 

background  for  energies  greater  than  -  20  eV,  but  highly  sensitive  for  lower 

energies.  Energy  lost  to  background  electrons  increases  with  increasing 

numbers  of  these  electrons,  however,  V  remains  nearly  constant  until  the 

_2 

fraction  approaches  10 

Finally,  we  note  that  the  results  of  this  study  are  sensitive  to  the 
magnitude  of  the  K-shell  cross  section,  when  included.  Unfortunately,  no 
experimental  data  on  K-shell  cross  sections  or  V  (for  0+)  exists  so  the  choice 
is  not  unique. 


1.  S.  Slinker  and  A.V.  Ali ,  Proceedings  of  the  SDIO/DARPA  Annual  Propagation 
Review,  Alburquerque,  New  Mexico  (1986). 

2.  S.  Slinker  and  A.V.  Ali,  "Electron  Energy  Deposition  in  Atomic  Oxygen", 

NRL  Memorandum  Report  5909,  Washington,  D.C.  (1986). 

3.  R.D.  Taylor,  S.P.  Slinker,  and  A.V.  Ali,  "Electron  Energy  Deposition  in 
Atomic  Nitrogen",  NRL  Memorandum  Report  6087 ,  Washington,  D.C.  (1987). 

A.  R.D.  Taylor,  S.  Slinker,  and  A.V.  Ali,  Proceedings  of  the 

SDIO/DARPA/Services  Annual  Propagation  Review,  Monterey,  California 
(1987). 

5.  S.  Slinker,  R.  Taylor,  and  A.V.  Ali,  Bull.  Am.  Phys.  Soc.  32,  1947  (1987). 

6.  A.V.  Ali,  S.  Slinker  and  R.D.  Taylor,  BEAMS  ' 88  Seventh  International 

Conference  on  High-Power  Particle  Beams,  Karlsruhe,  Germany  (1988). 

7.  S.P.  Slinker,  R.D.  Taylor,  and  A.V.  Ali,  J.  Appl.  Phys.  63,  1  (1988). 

8.  R.D.  Taylor,  S.P.  Slinker,  and  A.V.  Ali,  J.  Appl.  Phys.  64,  982  £1988). 

9.  R.D.  Taylor,  A.V.  Ali,  and  S.P.  Slinker,  "Energy  Deposition  in  0  by  High- 

Energy  Electron  Beams",  J.  Appl.  Phys.  (in  press). 

10.  H.A.  Bethe,  Handbuch  der  Physik,  Vol.  24,  p.  273,  Springer,  Berlin  (1933). 

11.  R.D.  Taylor  and  A.V.  Ali ,  "Excitation  and  Ionization  Cross  Sections  for 
Electron-Beam  Energy  Deposition  in  High  Temperature  Air",  NRL  Memorandum 
Report  6013 ,  Washington,  D.C.  (1987). 


24 


Transverse  Two  Stream  Instability* 

Glenn  Joyce  and  Martin  Lampe 

Beam  Physics  Branch,  Plasma  Physics  Division 
Naval  Research  Laboratory,  Washington,  D.C. 

INTRODUCTION: 

A  relativistic  electron  beam  propagating  through  a  diffuse  plasma  can 
interact  strongly  with  the  plasma  resulting  in  a  change  of  the  properties 
of  the  beam.  The  beam  can  generate  electrostatic  Wakefields  in  the  plasma 
causing  the  plasma  to  act  as  a  transformer  and  move  energy  from  one  region 
of  the  beam  to  another.  This  aspect  of  the  wakefield  effect  is  primarily 
axisymmetric  and  has  been  treated  extensively  in  the  past.  Another  aspect 
is  the  interaction  of  nonaxisymmetric  perturbation  of  the  plasma  with  the 
beam.  The  result  of  this  type  of  perturbation  is  an  m  =  1  instability. 

For  a  beam  propagating  in  a  diffuse  plasma  of  constant  density,  the 
instability  is  absolute  in  the  beam  frame  and  has  a  large  growth  rate. 
However,  if  the  beam  is  propagating  in  an  IFR  channel,  the  instability  is 
convective  and  under  most  circumstances  the  perturbations  will  pass  through 
the  beam  before  the  disturbances  have  grown  large  enough  to  disrupt  the 
beam.  The  reason  for  the  change  from  absolute  to  convective  is  that  for  a 
channel  of  varying  density,  the  beam  electrons  oscillate  with  different 
frequencies  depending  on  their  radial  positions  giving  rise  to  phase-mix 
damping.  We  have  studied  the  transverse  instability  both  analytically  and 
numerically  using  the  ELBA  three-dimensional  particle  simulation  code. 


DISPERSION  RELATION: 

We  begin  with  the  linearized  cold  fluid  equations  in  the  electrostatic- 
magnetostatic  approximation  and  consider  m  =  1  perturbations. 

**nl  -»  -» 

sr +  v^o +  Vvi =  °’ 
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=  -4nq  (n  +  n^). 
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We  skip  the  details  of  the  derivation  and  write  the  result.  If  we 
consider  perturbations  of  the  type  ei(o>C  +  2z),  the  dispersion  relation 
becomes 
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average  over  betatron  frequency. 
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For  the  case  in  which  there  is  no  channel  w  .  =  »  ,  so  all  beam  electrons 

ch  p 

2  2 

have  the  same  betatron  wavelength.  Using  a^n^  =  an  nQ,  the  dispersion 
relation  becomes  (v  -*  0), 

2 

t£_  -  1  1 

2  ~  1  ~  2  2  2’ 

(/  Z  1  -  2V2‘ 

P  3 

which  is  the  result  of  Yu  and  Sharp  (private  communication). 

For  the  case  of  a  beam  in  an  IFR  channel  and  a  diffuse  plasma,  the 

channel  density  varies  with  r  leading  to  a  spread  in  betatron  wavelength. 

Also  n  .  »  n  ,  so 
cn  o 

abnb  -  anno  *  ach"ch  ■  an”o/(1-£ch>- 


We  write  Sr  =  2/2^  and  h  =  2^/2^.  Then 


* 1 1 1  *  C°h)  ‘-lb  Cs2) 


which  makes  the  dispersion  relation  look  like  a  hose  dispersion  relation. 
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If  ve  take  the  spread  mass  model  for  a  Bennett-like  channel,  ve  can 
estimate  the  effect  of  the  instability  as  w^/io  =  (nQ/nch)  ^n/8.  Writing 

the  beam  length  Cmay  =  N  2nc/wQ,  where  N  is  the  number  of  plasma 

wavelengths  in  the  beam,  we  find 

(n  /nch)/(l-f)3n2/4  N 

Amplitude  =  e 

_3 

For  Delphi  N  =  28  and  n0/'ncj1  =  10  so  we  expect  the  amplitude  growth  to  be 
less  than  the  beam  radius. 

PARTICLE  SIMULATION: 

Ve  have  simulated  the  instability  with  the  ELBA  code  considering  beam 
propagating  in  diffuse  plasmas  with  and  without  IFR  channels.  For  these 
runs,  the  simulation  parameters  were  1=1  kA,  y  =  4,  and  a^  =  1  cm.  The 

simulations  were  for  either  a  constant  density  background  plasma  or  a 

channel  with  a^  =  1  cm,  and  fch  "  -5  and  a  constant  density  background 
plasma.  The  beam  and  plasma  densities  are  chosen  such  that  =  104  cm  and 
Ap^  =  37  cm. 

The  simulation  results  for  the  constant  density  case  are  shown  in  Figs. 
1-4.  For  a  constant  background  plasma  after  some  propagation  distance,  a 
purely  growing  mode  appears.  Figure  1  is  a  plot  of  the  centroid 
displacement  as  a  function  of  z  for  constant  £.  From  this  we  have  measured 
a  frequency  of  2  =  .91.  The  growth  rate  can  be  determined  from  Fig.  2 
which  is  a  plot  of  centroid  displacement  as  a  function  of  C  for  constant  z. 
The  growth  rate  determined  from  this  plot  is  =  1.6.  According  to  the 
dispersion  relation,  if  2  =  .91,  =  1.4. 

For  the  Bennett  channel,  the  growth  is  much  less  violent  and  appears  to 
stop  before  the  beam  oscillations  become  as  large  as  the  channel.  Figure  3 
is  similar  to  Fig.  2,  and  Fig.  4  is  similar  to  Fig.  2  except  in  this 
simulation,  a  channel  of  strength  .5  was  included.  The  effect  of  the 
channel  is  clearly  to  limit  the  growth  of  the  instability. 
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CONCLUSIONS: 

A  relativistic  beam  propagating  through  a  diffuse  plasma  is  unstable  to 
transverse  perturbations.  For  a  plasma  of  constant  density,  the 
instability  is  quite  severe  and  can  rapidly  disrupt  the  beam.  The 
inclusion  of  an  IFR  channel  with  varying  density  can  change  the  nature  of 
the  instability  and  cause  the  instability  growth  to  stop  before  while  the 
beam  is  still  in  the  channel.  Ve  have  presented  and  analytic  theory  and 
supporting  simulations  of  this  effect. 


♦Work  supported  by  the  Defense  Advanced  Research  Projects  Agency,  ARPA 
Order  No.  4395,  Amendment  80,  monitored  by  Naval  Surface  Warfare  Center. 
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